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Abstract 

We report the persistent current in two coupled rings which form a character "8" structure. 
We showed that the magnetic fluxes modify the global phase of the electronic wave function for 
multiple connected geometry formed by the coupled rings. We obtained an exact solution for the 
persistent current and investigated the exact solution numerically. For two large coupled rings 
with equal fluxes, we found that the persistent current in the two coupled rings is in fact equal to 
that in a single ring. For opposite fluxes the energy has a chaotic structure. For both cases the 
periodicity is h/e. This results are obtained within an extension of Dirac's second class constraints. 
This theory can be tested in the ballistic regime. 
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I Introduction 

In Quantum Mechanics, the wave function is path dependent and is sensitive to the 
presence of a vector potential caused by an external magnetic flux. In a closed geometrical 
structure, such as a ring the wave function is changed by a measurable phase [1], causing 
all the physical properties to become periodic [2]. When a mesoscopic ring of normal metal 
is pierced by a magnetic flux $ [2], the boundary conditions are modified, leading to a 
famous theorem of periodic properties with the flux period $o = h/e and to a remarkable 
phenomena [3] of a non-dissipative persistent current [3-7]. 

One way to classify the closed geometrical structure is by using the number of holes 
formed on the closed geometrical structure. For a closed surface, the number of holes 
formed thereon is often referred to as a genus number g [9, 11]. For example, a genus 
number g = \ describes an Aharonov-Bohm ring geometry, while a genus number g = 2 
describes two rings perfectly glued at one point to form a character "8" structure. 

In this paper, we report the first exact solution for multiple connected geometries, such 
as a geometry with two holes two rings perfectly glued at one point to form a chcirac- 
ter "8" structure, or connected rings). The geometry modifies the global properties of the 
wave function, and the presence of magnetic fiuxes generates persistent currents with com- 
plicated periods. We present an exact analytical solution for the eigenvalues and compute 
the persistent current for two coupled rings with a character "8" structure for two different 
fiuxes. We solve the problem by modeling the gluing of the two rings using Fermionic 
constraints with anti-commuting Lagrange multipliers, which can be viewed as a resonant 
impurity strongly coupled to the two rings . 

The analjd;ical results are investigated numerically. When the two fiuxes on both rings are 
the same, we find a simple relation between the single ring [g = 1) current, I^^^^\flux] N), 
and the double ring {g = 2) current, I^^^'^\flux; N). At T = 0, we define I^^^'^\flux; N) = 
r{N)I^^^^\flux; N), where r{N) is a ratio between the two currents. The ratio r{N) is a 
function of the number of sites N and obeys r{N) — > 1 for A?" ^ oo. 

The plan of this paper is as foUowings in chapter II we present the exact analytical 
results for the two rings perfectly glued at one point to form a character "8" structure. In 
chapter III we present the numerical results for the two coupled rings for equal and opposite 
fines. Chapter IV is devoted to the discussions: We show that the method presented can 
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be extended to two rings with a finite with. 

Using the one dimensional Persistent current (single channel) given in ref.9 in the presence 
of a 2Kp impurity scattering given in ref.9 we compute the current to the multichannel case. 

We show that the theory developed for the Ballistic regime might be able to explain the 
experimental results for the GaAs/GaAlAs coupled rings [13]. 

We also present an generalization of the Dirac's constraint method for a genus g > 2 
multi coupled rings. 

II Exact Solution for two rings perfectly glued at one point to form a chciracter 
"8" structure 

Dirac has shown [14 ] that for the second class constraints the Poisson brackets are 
replaced by the Dirac [14] brackets. For an even number of constraints q^. with non zero 
Poisson brackets., the equations of motions are governed by the Dirac [14] brackets which 
replace the Poisson bracket {A,B} by ,{A,B}d — {^,B} — Ylirr'{-^^^r}crr'{(lr',B} . The 
matrix €„' is given in terms of the constraints qr , ^ Crr'{qr',Qs} — Sr,s- To obtain a Bosonic 

r' 

theory one replaces the Poisson bracket by the commutators, {, } = ]• 

We propose that for the second class Fermionic constraints the following modification. 
Given two Ferminonic constraints Q, which obey non-zero anticommutation relations, 
\_Qi + = QQ"^ + Q^Q 7^ 0. We find that the Dirac bracket is replaced for a Fermionic 
operator Op and the hamiltonian Hq by: 



0,Ho 



D 



0,Ho 



0,Q 



( [Ql Q])'' [Q: Ho] - [O, g] ^ ( [Qt, Q] [g, Ho] 

Therefore the new Heisenberg equation for any fermionic operator Op will be given by: 

Op,Ho . For the remaining part we present the derivation and applications of 
J D 

this new result. 

We consider the Hamiltonian Ho for two spinless Fermionic rings in the absence of a 
magnetic flux. The rings obey periodic boundary conditions. For each ring, the point x is 
identifled with the point x + L. The two coupled rings with the character "8" structure (i.e. 
g — 2) are obtained by identifying the middle point x — L/2oi the flrst ring with point x — Q 
of the second ring, i.e. Ci(L/2) = 6*2(0) and C^{L/2) = C^(0). This identification is equiv- 
alent to two Fermionic constraints, Q = Ci{L/2) — 6*2(0) and = C^(L/2) — 6*2^(0). 
Since the constraints are Fermionic, they can be enforced by using anti-commuting La- 
grange multipliers, /x"*" and /x. Following ref. [14], we introduce the Hamiltonian with the 
constraints, — H + /I'^Q + Q'^/i. The unusual physical meaning of the anti-commuting 
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Lagrange multipliers can be viewed as a Fermionic impurity [14], which mediates the hop- 
ping of the electrons between the two rings. This method is simpler in comparison with 
the method based on matching boundary conditions for the wave function explained in the 
discussions paragraph. The two rings of length L are threaded by a magnetic flux where 
a = 1,2 (for each ring). In order to observe the changes of the constraints in the pres- 
ence of the external flux, we perform the following steps. In the absence of the external 
flux $0,, the annihilation and creation fermion operators obey periodic boundary conditions 
Ca{x) = Ca{x + L) and C'^{x) = C^{x + L), where a = 1,2. The genus g = 2 is imple- 
mented by the Fermionic constraints Q = Ci{L/2) - ^2(0) and Q+ = Ct{L/2) - C^{0), 
and Hamiltonian Hq = — t X]a=i Si—V^^^i^Q (^)^a(^ + a) + h.c.]. The length of each ring 
is L = Nga, where A^^ is the number of sites and a is the lattice spacing. When the 
external magnetic flux is appled the Hamiltonian Ho is replaced by H. The Hamilto- 
nian H is obtained by the transformation Ca{x) exp[i^ A{x'; a)dx']Ca{x) = ^ipaix) and 
^ai^) ~* 6xp[— 2-^ v4(a;'; a;)(ix'] = ip'^{x). Here A{x]a) is the tangential compo- 

nent of vector potential on each ring. The relation between the flux and the vector potential 
on each ring is A{x; a)dx = (fa- 

The flux on each ring a = 1, 2 gives rise to a change in the boundary conditions, 
ipa{x + Nsa) = V^„(x)e^'^" and ^p+{x + A^,a) = V^+(x)e-*'^% where cpa = 27r(^) = 2tt^ = 
2TT0a- This boundary condition gives rise to a normal mode expression for each ring, tpa{x) = 
S^lio^ ^ similar expression for ip^{x). The "momentum" K{n,ipa) 

is given by, K{n, ip^) = -^{n + 0a) where n = 0, 1, . . . , — 1 are integers with N = Ng, 
and (pa = 2'K0a- In the momentum space, the Fermionic operators ipa{n) and ip^lm) obey 
anti-commutation relations, [ilJa{n),ip^{m)]^ = Sa,i3Sn,m- The Hamiltonian for the two rings 
in the presence the external flux takes the form, 

{Ns-l)a Ns-1 

H = -tJ2 XI ^«(^)^"^^ + ") + ^-^- = 5Z 5Z ^(^''^°)^-^'^)^-^'^) 

a=l,2 x=0 a=l,2 n=0 

where e{n,ifa) = — 2tcos[^(n + 0a)] are the eigenvalues for each ring. The Hamil- 
tonian in eq. [1] has to be solved together with the transformed constraints ,Q = 

The wave function for the genus g = 2 problem is given by the eigenstate \x) of the 
Hamiltonian in eq. [H which in addition satisfies the equations Q\x) = and Q^\x) = 0. 
The constraint conditions are implemented with the help of the anti — commuting Lagrange 
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multipliers fi and /i"*". The Hamiltonian Ht with the constraints takes the form, 
HT = H + fi+Q + Q+fi 

The Lagrange multipher are determined by the condition that the constraints are satisfied 
at any time. Therefore, the time derivative satisfies the equation, Q\x) = Q^\x) = at 
any time. We will use the notations, [A, S]+ = AB + BA and [A, B] = AB - BA. The 
Heisenberg equation of motion for the constraint Q is, 

ihQ = [Q, Ht] = [Q, H] + [Q, /x+Q + Q+fj] 

= [Q, H] + [Q, /i+]+g - fi+[Q, g]+ + [Q, Q+]+/i - Q+[Q, 

= [Q,H] + [Q,Q+Ui, (2) 

The rest of the anti-commutators in eq. [2] vanishes. The anti-commuting Lagrange mul- 
tipliers satisfy, [Q,fi'^]^ = [Q,fi]^ = = [Q^,^]^ = 0. Since the constraints are 
fermionic, we obtain that they obey [Q,Q^]-^ = [Q^,Q]+ = 2. Therefore, the constraints 
are second class constraints [14]. From the condition Q\x) =0 and eq. [21 we determine the 
Lagrange multiplier field /i, 

fi = -[Q+,Qn'[Q,H] = -l[Q,H]. 

The field /i"*" is obtained from the equation Q~^\x) = 0, 

^+=[Q^Q+]-l[Q+^H] = l[Q+,H]. 

The Hamiltonian Ht with the constraints and the Lagrange multipliers are used to 
compute the Heisenberg equation of motion for any Fermionic operator, O. (The Lagrange 
multipliers anti-commute with any Fermionic operator, i.e. [0,/i]+ = [0,p'^]^ = 0.) 

th^ = [0, Ht] = [o, H] + [0, /x+g] + [6, g+]^ 

= [6, H] + [6, fx+]+Q - /i+[6, g]+ + [6, g+]+/i - g+[d, 

= [d,H]-[d,Q]+^^^-[d,Q+]^i (3) 

We substitute in eq. 3 the solutions for the Lagrange multiplier fields and obtain a new 
equation of motion with a new commutator, which resemble the classical Dirac brackets 
[14]. 

th^ = [d,HT] = [d,H] - [6,g+]+([g+,g]+)-i[g,/7] - [d,Q]^i[Q,Q+]^)-'[Q+,H] 

= [0,H]n (4) 
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Eq. m shows that the Heisenberg equation of motion is governed by a new commutator, 
[0,H]d- The equations Q\x) = and Q'^lx) = are inconsistent with [Q,Q^]\x) 7^ 0. 
The new commutator resolves the inconsistency problem, [Q, I x) = 0! We will use 
this new commutator to compute the Heisenberg equations of motion for the creation and 
annihilation Fermionic operators ipa{x,t) and ip'^{x,t), where a = 1,2. 

ihMx) = [Mx),h]d = [Mx),h]-^[Mx),q+]4q,h] 

= -t[ilJo,{x + a) + tpo,{x - a)] - ^[6a,iS^^L/2e"^^ - 5a,24,o] 

i-t){e-'^nM^ + a)+M^- a)] + e-^-^[u\ + «) + u\ " «)]} (5) 

The ground state wave function is obtained from the one electron state, \x >= 
^a=i2^x^o~^^^ ^aix)^lj^{x)\0 >, giveu in terms of the site amplitudes Za{x). Using a 
normal mode momentum expansion, fa{n), i.e. Za{x) = -^J2n=o^^^^"'''^°'^^fcii^)^ 
the following equations for the eigenvalues A and the amplitudes in the momentum space 
fa{n), 

(A - e{^ + = E <^ + - ^ E <^ + ^2)/2(^) (6) 

n=0 n=0 

and 

(A - 6(£ + = ^ E <^ + ^2)/2H + 1^ E <^ + (7) 

Tl=0 71=0 

We diagonalize eqs. [6] and [7] by linear transformations, S'i((^i,A) = —Yl!t=o + 
and S'2(<^25A) = —Yld=^^{^^'^'i)f'i{^)- ^ result, we obtain the equation, 

= 0, where the matrix M is given by M = ^ x • Here, 

we define aL^^((^„,A) = Ai''''^"''(<^„, A) + Ai°'^'^''(<^a, A) and Ai"^((^a,A) = aL'''"'"''(<^q, A) - 
Aa'^'^\0a, with the even and odd representations given by, Aa^^"'\0a, ^) = 

1 v^(jV-2)/2 e{2m+^c) i AHd),',^ \\ _ 1 v^(jV-2)/2 e{2m+l+0c) ^ pomniltp 

2ArZ^m=0 X-e{2m+M IV^a, - 27V Z^m=0 A-.(2m+l+<^<,) " "'^ COmpUte 

det M = and obtain the characteristic polynomial which is used to determine the 
eigenvalues A. 

2[A[--\^„ X)Ai"''\^2, A) + A[°''\^„ X)Ai--\^„ A)] + [a[^\^„ A) - a(+)(^2, A)] = 1 

(8) 




Eq. [8] is our main result for the genus g = 2 case. We observe that the matrix M is 
symmetric and the eigenvalues are real when the fluxes are equal, i.e. 0i = 02, or opposite, 
i.e. ifi = —02- For other cases, the eigenvalues can have imaginary parts, thereby giving 
rise to non-conducting states. 

Ill Numerical Solution 

We have numerically solved the secular equation [HI To compute the current, we 
sum over the current carried by each eigenvalue X{0i,02) using the grand-canonical 
ensemble. The current in each ring a = 1,2 is given by, Ia~^\(pi, ^2) = 

- 4: [^(^^- ^^)]^( ^'^t::L!:r'^ ) where is the Fermi Dirac 

function with the chemical potential Efermi and temperature T. The current is sensitive 
to the number of electrons being either even or odd. We use the grand-canonical ensemble 
and limit ourselves to a situation with even numbers of sites and a zero chemical poten- 
tial, i.e. Efermi = (which corresponds to the half-filled case). In order to have a perfect 
particle-hole symmetry, we will restrict the analysis to the special series for the number 
of sites being A^^ = 2, 6, 10, 14, 18, ... , 2m + 2, where m = 0, 1, 2, 3 . . .. For this case, we 
find that, when the fluxes are the same in both rings, the current for g = 2 has the same 
periodicity as that of a single ring, i.e. 

/(9=2)(^$ + <|,p) = /(9=2)($). At temperatures T < 0.02 
Kelvin, the line shape of the current as a function of the flux is of a sawtooth form (see 
figure 16). For other series A^^ 7^ 2m + 2, the periodicity of the current is complicated. 
Using the experimental values given in the experiment [11], we estimate that the number 
of sites in our model should be in the range of A^s = 50 ~ 150, the hopping constant 
should he t = t. — — 7 ~ 0.01 eV and the temperature in the experiment should be 
T = 0.02 Kelvin. Using these units, we obtain that the persistent current is given in terms 
of a dimensionless current, I (see figure 16 and figure Ic) with the actual current value, 
j(3=2) ^ J X Q 92 X 10-4 Ampere. 

a) Equal fluxes for g=2 

For this case the secular equation is simplified and takes the form of 
4[A(^^'^")(^,A)A(°'^'^)(^,A)] = 1. 

For A^s = 2, we solve analytically the secular equation. We find that the eigenvalues are 
given by \{n,Lp]N = 2) = r{N = 2)e{n,yD;N = 2), where e{n,yD,N = 2) = -2tcos[^(n + 
(p)], and n = 0,1 are the single ring eigenvalues. The value for r{N = 2) is r{N = 2) = 



7 



To find the eigenvalues for otlier number of sites, A^^ = 6, 10, 14, 18, 22, 26, 30, we numerically 
find the relation, A(n, yj; A^) = r{N)e{n, ip; N), where = 0, 1, . . . , — 1 and e(n, ip; N) = 
are the single ring eigenvalues. The function r{N) is given in figure la. 
This figure shows that the function r{N) reaches one for large A^. Using the function r{N) 
given in figure la, we compute the current for the g = 2 case as a function of temperature, 

/(-^)(^; N-T) = - IHN) . e{n, ^- N)]F ^^''^tZZl'r"'' ^ 

Figure Ih represents the current for A^^ = 30 sites at two temperatures T = 0.02 and 
T = 20 Kelvin. In this figure, the current is given in dimensionless units / plotted as a 
function of the dimensionless flux f = (pa = [—0.5,0.5] {ipa = 27iipa = [—^j^r]). The solid 
line represents the single ring current and the dashed line represents the current for the 
genus g = 2 case. In figure 16, the ratio of the currents at T = 0.02 Kelvin is r{N = 
30, T = 0.02) = 0.979. 

Figure Ic shows that the currents at T = 20 Kelvin, are in the range of 7 nA and the 
reduction of the current is larger in comparison with the T = 0.02 Kelvin case given in figure 
16. 

b ) Two coupled rings with opposite fluxes , i. e. (pi = —02 

For Ng = 2, the eigenvalues are the same as the one obtained for the same flux 
case. For A^^ = 6, 10, 14, ... , 2m + 2, we solve the secular equation given in eq. [8] 
and compute the eigenvalues. In figure 2a, we plot the total energy as a function of 
the opposite fluxes at T = 0.02K for 30 sites, E^s=^\-0,0, N, = 30,T = 0.02, iT) = 
Yln=o~^['^i~'fi' Kjfl '^■^y""'' )] . The total energy dependence on the opposite flux 

is chaotic due to the interference between paths which encircle zero and non zero fluxes, 
caused to the common point between the rings which acts as an impurity. In addition we 
observe periodic oscillation with the fundamental period $o (see figure 2a and 2b). For 
comparison, we show in figure 2b the total energy for equal fluxes which is parabolic and 
the current is linear (for small fluxes). 

The effect of a finite width in the ballistic regime will give rise to a multichannel Persistent 
current which can be analyzed using the discussions for the multichannel case presented in 
the next chapter. Therefore the energy will be given by a sum of energies with respect 
the different channels. Therefore we expect that the multichannel effect will give rise to a 
smooth function of the energy as a function of opposite flux . 
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IV Discussion 

a) Possible experimental application of our results. 

At this stage it is not clear if a genus g—2 experiments exists in the ballistic regime where 
our theory can be applied. The closest experiment which might be relevant to our theory 
is the experiment presented in ref.l3. The author of ref. 13 make the statement (see the 
last sentence in their paper) that from the theoretical side a model for the ballistic regime 
is needed for a direct comparison with the experiment. They make this statement based on 
their results which indicate that the ratio between the single ring current and the current 
in 16 coupled rings is close to one! The authors of ref. 13 say that the conditions in their 
experiment are not in a clear diffusive regions therefore we can try to apply our ballistic 
theory to explain the experiments. 

The results presented in ref. [13] are based on a system of 16 GaAs/GaAlAs coupled 
rings. At this stage we have only results for two rings: We find for two rings the ratio 
I^^^'^y I single-ring — 0.987. Sincc the experiment was performed on 16 rings we use a scaling 
argument in order to extrapolate the results to 16 rings. For two rings plus a scaling 
argument we obtain, r = he-rings / Lingie-Hng ~ [1''^^'^^ Lingie-ring]'^ = [r{T = 0.02, Ns = 
50)]^ = [0.987]"^ = 0.95. The value r = 0.95, is in the range of the experimental observation 
reported [13]. In the experiment the rings are connected trough arms of length of the order 
of the wavelength which can be approximated by point contact between the rings. 

Next we evaluate the amplitude of the current. Since at low temperature and large 
number of sites we have found that the ratio is close to one we will consider the amplitude 
for a single ring. Using the experimental values for the Fermi velocity vp — 3.16 • 10^ • ms"^, 
Fermi wavelength Xp = 3.5 • 10~^m and the ring perimeter L = 1.2 • 10~^m we compute the 
persistent current for a single ring. Using a model of a single conducting channel we find at 
that the amplitude current at T = is given by, ^ = A.2nA. Using the velocity and Fermi 
wavelength we obtain the effective hopping constant used in our simulations, t — 2asin^Kj-a) ' 
using Kp — 1^ with the effective lattice constant a — Xp/^. This gives us that the effective 
length used in the simulation for 30 sites was Lgimuiation — 30 • Xp/A: — 2.6 • 10~^m. Our 
simulation show that the current for 30 sites at T = 0.04 K was / — 10~^ in dimensionless 
units which corresponds to a current I ^ I ■ 0.92 ■ lO^^A. Using the value of the hopping 
constant t = 1.9Joule = 0.0126.^. we find for Ng a current , I{Ns = 30) = = ir=m ' 
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10-6 = 92nA 

In figure 3 we show the current dependence on the length of the ring for single and double 
rings at different temperature. At T=0.02 K the current decreases hnearly with the length. 

This allows to use at T=0 the linear relation between I{L = 1.2 ■ lO^^m) = A.2nA and 
HLsimuiation = 2.6 • lO'^m) = 92nA which obey the relation, ^^^"7"'°*'°"^ = j — 



-'simulation 



b) The effect of disorder on the Persistent current. 

We explain the discrepancy of a factor of 10 between the experimental result and the 
theoretical calculation using the 2/ri?impurity scattering for a multichannel one dimensional 
rings at zero temperature. 

Contrarily to the transport current the Persistent current is without dissipation. The 
persistent current is determined by the time derivative of the zero mode coordinate (the 
macroscopic phase of the wave function). In ref.9 (Schmeltzer and Berkovits) it was shown 
that the effect of a 2Kp impurity scattering gives rise to an enhancement of the kinetic mass 
and as a result the Persistent current is suppressed. Using the results given in figure 1 (ref.9 
Schmeltzer and Berkovits) we find that the suppression factor for one channel was 1 — bW^, 
where b — 38.29 and Wis the 2Kp impurity scattering. The value of W can be obtained from 
the transport value 1^, h — ^- lO^^m. Prom the Boltzman equation we will extract W, = 
< (1 — Cos{9)) >. The angular average is determined by the ratio of the width 
and length of the ring. We have, < 1 - Cos{9) >= J^I^{1 - Cos(^)) f + - Cos(^)) f 

with 9o = where d = lO^^m is the width of the ring (which is actually a square with 
the length Le// — . Solving this equation we find that scattering potential is given by, 
T = (7fli)^/' = 0.16. 

Using the result given in ref.9 (Schmeltzer and Berkovits) we compute the Persistent 
current for a single ring (with one channel) in the presence of the 2Kf potential f = 0.16. 
We find ^{1-bW'^) = 4.2nA-(l-38.26-(0.16)2) = 0.09nA which is smaller than the current 
observed in the experiment, 0.4nA. This discrepancy suggest the possibility of a multichannel 
effect. For a single ring the wave function is given by, 'i>{x,y) — Yln=T '^n{x)Xn{y) where 
Xn{y) — (i)^^^'5'^'T'('T'f ) is the standing wave in the transversal direction of the ring < y < d. 
In the absence of disorder the problem is replaced by n independent one dimensional channels 
with a shifted energy ^{n^)'^ and operators i/Jni^) ,'ijjl^{x). As a result the current will be 
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given by, I(nodisorder) = ^YZ=Ti^ - i'^^^TY^'^- The effect of the 2Kp impurity 
scattering will replace the current by, 

J = 2^ ^^—(1 - (nf ^)2)V2[i _ i,w'^l - (nf If ( we have replaced b by 6(1 - 

{^2^)^)^^ , in agreement with ref.9 where it was shown that 6 is a function of the Fermi 
velocity for each channel). Using the extracted values ofb and W we find from the condition 
[1 — bW'^{l — ^ that the maximum number of conducting channels is given by, 

nmax — 8. Using the value of Umax — 8 we compute the current and find that the Persistent 
current is given by I — OAnA which is in agreement with the experiment in ref.13. 

the constraints for two rings which have a finite width. 

For this case we have to do the gluing for a ring of width d. The transversal direction is y 
and xis the direction of the one dimensional cchannel. The constraints used previously are 
modified in the following way: Q{x — Q,y — d + e) — Ci{L/2, y — d + e') — 6*2(0, y — d + e") 
and Q+(0, y = d + e) = C+(L/2, y = d + e') - Ct{Q, y^d + e"). 

Where e' + e" — larm and larm — d represents the length of the arm which connects the 
two rings. As a result the common part between the rings will be given by the connecting 

arm larm — d. 

The second class constraints is given by, 

[Q{x ^Q,y^d + e'), Q+{x' ^Q,y = d + e")]+ = 2d{e' + e" - larm) 

Using the second class constraint we obtain the equation of motion for the two rings spinor 
Ca=i,2{x,y) 

ih'^^ = [C^{x,y),HT\D = [Ca{x,y),H] 

-f dzf dz'[Ca{x, y), g+(0, y^d + ^]+([Q+(0, y^d + z), Q(0, y = d + z']+)-'[Q{0, y = 
d/2 + z'),H] 

-f dzf dz'[C^{x,y),Q{0,y^d + z] + {[Q{0,y^d + z),Q+{0,y^d + z')U)-'[Q+{0,y^ 
d + z'),H] 

In the absence of disorder we can simplify the expression for the constraints if we introduce 
a contact point y — do chosen such that,d < do < larm)- As a result we obtain a constraints 
for each channel, Qnix = 0) and Qmix = 0). We have the conditions for the second class 
constraints, [Qn, Qm]+ = '^S{n,m)- As a result the current for two rings will be similar as the 
one dimensional case . The difference between the current for different channels in two rings 
will be determined by the transversal shift potential energy, ^{n^)^. We have repeated our 
simulation for different transversal potentials with n < Umax — 8 (the number of propagating 
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zero modes ) and find that the ratio between the double and single rings remains one. 



c) Generalization to many coupled rings 

In order be able to study a system of many coupled rings we have to find a method 
which can be applied to many coupled rings . It seems that it might be preferable to replace 
the constraint of the single particle operator with the constraints defined in terms of the 
electronic densities and currents. This method can be combined with the basic technique of 
matching the boundary conditions used for for solving Quantum wires problems. 

At the common point of the two rings at x=0 the constraint should gives rise to equal 
densities ,pi (x = 0) = p2 (x = 0). Formally this condition is enforced with the help of the 
scalar field a^d {x) which plays the role of the Lagrange multiplier. For g-coupled rings 
we introduce a set of scalar potentials Cq^'^'*..., ..Cq^"^'^"^^ for which a statistical annealed 
average has to be performed. 

The hamiltonian for "g^" coupled rings with p„ {x) electronic density and </?i </?g fluxes , 

= H i^,...^,) + E^„:5 4"'"^'^ (Pn (x - (n + 1) L) - {x - nL)) 
To demonstrate this method we repeat the case g—2 two ring case. 
For a symmetric configuration with the common point at x = 0, we fold the space of the 

first ring from [— i^, 0] to [L, 0] such that the space of the two ring is restricted to < x < L. 

The wave function for the coupled rings Ze{x) with eigenvalue E is given as a spinor with 

two components Zt (x), t — 1,2. The Schrodinger equation for the two rings in the presence 

of the field aoS (x) is, 

V 



27r 

-dx - i-rfi ) + «o(^ (x) 



Zi (x) = K^Zi (x) 



.27r 



- [dx- ^^V2 ) - ao5 (x) 



Z2 (x) = K^Z2 (x) 



Where K = with the energy E. The current / [(pi, (p2; Oq] is a function of cpi, (p2 and 

The eigenfunctions for this problem are obtained by matching the boundary conditions: 
a)the continuity of the wave functions at the boundary x—0 , Zi{x) — Zi{x + L) and 
Z2 (x) = (x + L) 

b) The discontinuity of the derivative at x—0 for each ring, 
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{-d^ - ifcpi) Zi {x = -e) - {-d^ - ifcpi) Zi (x = £) = aoZi {x = 0) 

{d^ - ifip2) Z2 (X = -e) - {-d^ - i^ip2) Z2 (X = £) = -00^2 (x = 0) 

(The change of sign of the flux in ring one a = 1 is a result of folding the space from 

[-L,0] to [L,0].) 

As a result the current / (/92; Oq] is a function of ipi, and the scalar potential Oq. 
The physical current will be obtained after averaging over the constraint Qq, / ((^1,(^2) = 
lim J^,^I{(pi,(p2;ao). The need for the additional average makes the solution for the 

A— »oo 

persistent current more involved. Therefore the resulting method is more complicated in 
comparison with Dirac's method used in the first part for two rings. From other hand for 
for many rings the new method seems to he preferable to the Dirac's method. This becomes 
clear once the electron-electron interactions and disordered are considered, this new method 
allows the use of field theory methods such as BOSONIZATION and Renormalization Group 
[9,10]. 

Summciry 

In this paper, we have introduced a method which solves the problem of the global phase 
of the wave function for geometrical structures with holes, i.e. high genus materials. This 
method is applicable to a variety of mesoscopic systems where coherency of wave function 
is important. 

We have found an exact solution for the persistent current in two coupled rings. By 
numerical calculations, we have computed the current dependence on the flux, temperature, 
and the number of sites. This theory might be tested in coupled rings for equal and opposite 
flux in the Ballistic regime. A possible explanation of the experiment in ref.l3 has been 

proposed. 
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FIG. 1: (a) The ratio of the double to single ring currents I{g = 2;N)/I{g = I; N) = r{N); (b) 
The single ring (solid line) and the double ring (dashed line) currents for A'^g = 30 at T = 0.02 
Kelvin; and (c) The single ring (solid line) and the double ring (dashed line) currents for Ns = 30 
at T = 20.0 Kelvin. 
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FIG. 2: (a) The total energy for opposite fluxes, (/) = ipi = —ip2 for 30 sites at T = 0.02 Kelvin 
^(9=2)(_0 0. jv^ = 30,T = 0.02K) ;and (b) The total energy for equal fluxes ,0 = </?i = (p2 
^(9=2) (</), </); Ns = 30, T = 0.02ii:) 
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FIG. 3: The current derivative with respect the flux at (j) « 1 multiplied by the number of sites 
as a function of N at different temperatures-/' • for a single and a double ring 
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